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Abstract. Through polarization observations H 2 masers are excellent probes of magnetic fields in the maser region. Magnetic 
field strengths, such as those in the H2O masers regions of the envelopes of late-type stars and star-forming regions, are typically 
determined using a direct relation between the field strength and the observed circular polarization. Here it is shown that velocity 
and magnetic field gradients along the maser have a significant effect on the field strengths obtained from circular polarization 
observations. Due to velocity gradients the actual magnetic field strength could be up to 100% higher than the field strength 
derived from the observations. Additionally, when a magnetic field gradient is present, the resulting circular polarization derived 
is caused predominantly by the average magnetic field in the unsaturated maser core. Measurements of the fractional linear 
polarization are not affected by velocity or magnetic field strength gradients, though changes in the magnetic field angle along 
the maser do quench the linear polarization intensity when the maser saturates. 
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1. Introduction 

Astrophysical H2O masers, like other maser species such 
as OH and SiO, are excellent probes of magnetic fields 
in a variety of interesting regions. The magnetic field 
strength is derived from observations of maser polariza- 
tion caused by Zeeman splitting. Polarization observations 
of the 22 GHz H2O masers in the circumstellar envelopes 
(CSEs) around evolved stars have made important con- 
tributions to the understanding of the magnetic fields of 
evolved stars JvTemmings. Diamond. & van Langevelddiool 
IVlemmings et alJl2005l hereafter V02 and V05a). H?Q maser 
observations have also provided information on the strength 
and structure of the magnetic field in star - formi ng regions 
jFiebig & Guested ll989tlSarma et alJl200lL l2002t) and have 
provided uppe r limits for the f ield in the megamaser galaxy 
NGC 4258 llModiaz et all2005h . 

However, the exact relationship between the observed 
maser polarization and the magnetic field strength is not 
straightforward, leading to uncertainti es in the determ ined field 
strength. It was shown in iNedoluha & Watsod ( 1 19921 hereafter 
NW92) that the influence of the different hyperfine components 
of the H2O (6 16 - 523) rotational transition as well as of the 
degree of maser saturation is significant. The effect of maser 
saturation on the magnetic field strength determination can be 
as large as a factor of 2 (NW92). Additionally, velocity and 
magnetic field gradients might play an important role. While 
typical H2O maser spot sizes are between ~10 12 and 10 13 cm 



jReid&Moradll98lt llmai et aljf 19971) . the actual size of the 
maser region will be sever al times larger due to beaming (e.g. 
iGoldreich & Keelevlll972t) . As a result the maser path length 
could be as long as 10 14 - 10 15 cm. Across such a path length, 
velocity and magnetic field gradients likely influence t he ob- 
served polarization. In Vlemmings & van Langeveldel J2005L 
hereafter V05b) it was shown that velocity gradients alter the 
shape of the total intensity lin e profile. From this a nalysis and 
from earlier observations (e.g. lRichards et al J 19991) it is appar- 
ent that in CSEs velocity shifts of ~1 km s _1 are common. In 
this paper the effects of velocity and magnetic field gradients 
along the maser path on the magnetic field strength determined 
from circular polarization measurements and on the observed 
fractional linear polarization are examined. 

2. Background 

Similar to the analysis in V02 and V05b, the equations of state 
for the populations of the upper (616) and lower (523) rotational 
levels of the three strongest hyperfine components (F — F' = 
7-6, 6-5 and 5-4) of the 22 GHz H2O maser have been solved 
using the method described in NW92. The equation of state for 
the number density n(F, a, v) of the upper energy levels (F) is 



= A F (ys)-(T + T v )n(F,a,v) 

+R(F, F', a, b, v)(n(F', b, v) - n(F, a, v)) 
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Fig. 1. Normalized total intensity (I) and circular polarization (V) spectra for v t h = 0.6, 1.0 and 1.5 km s~'in the presence of a 
velocity gradient. From thick to thin the solid lines denote velocity shifts of AV m - 0, 0.5, 1 .0 and 1 .5 km s^'respectively except 
for = 0.6 km s _1 where AV,„ = 1 .5 km s~'is omitted. 



The equation for the number density n(F', b, v) of the lower 
levels (F') is similar but reversed. Here a and b denote the dif- 
ferent magnetic sub-states of each hyperfine component. The 
population levels are solved as a function of molecular velocity 
v and at different positions along the maser propagation direc- 
tion. The statistical weigh ts are designated by g F and g' F , and 
are from Kukulich (1969). The pump rate Ap(v s ) is assumed 
to be the same for the different hyperfine components and has 
a Maxwellian distribution. The results of the calculations de- 
pend only on the the ratio (Af - A' f )/Af, w hich is of the order 
of a few percent dAnderson & Watsonll993l) . Generally v s = v, 
but when velocity gradients are introduced later, the molecu- 
lar velocity of the maser pump v s changes along the maser 
path. The rate for stimulated emission R(F, F' , a, b, v) (here- 
after R) is calculated using the local maser intensity and the hy- 
perfine interaction coefficients as described in NW92. T is the 
decay rate for the molecular excitations. The cross-relaxation 
rate T v describes the reabsorbtion of previously emitted in- 
frared pump photons trapped in optically thick transitions of 
the maser system, which generate a newly excited molecule 
at random, within a Maxwellian velocity distribution {<p(v s )). 
Elastic collisions between the maser molecules and intermixed 
H2 molecules have not been included. The effects of these are 
described in detail in lElitzurl fl990) and are unimportant for the 
H2O masers discussed here. 



The maser intensity, which is solved using the radiative 
transfer equations from NW92, influences the level popula- 
tions through R. However, as it is itself dependent on the level 
populations determined with Eq. ^ the maser total intensity 
(I), linear polarization (Q) and circular polarization (V) are 
solved iteratively along the maser path. Assuming nearly one- 
dimensional maser propagation the beaming of the maser ra- 
diation is represented by the solid angle AQ. The emerging 
maser fluxes will thus be represented in 7bAQ, with 7b the 
brightness temperature. All calculations were performed with a 
brightness temperature for the radiation incident on the maser 
region of 7b AO = 0.1 K sr. It was verified in NW92 and V02 
that the results are insensitive to the chosen initial value. The 
calculations are performed for different thermal line widths 
(v t h) of the Maxwellian particle velocity distribution, where 
v t h ~0-5(T /100) 1 / 2 km s^with T the te mperature of the masing 
molecules in K. Nedoluha & Watsonl dl99ll) and NW92 have 
shown that when Eq.[0is solved without including the contri- 
bution of the cross-relaxation rate F v , the results for the line 
profiles, fractional polarizations and line width in the case of 
non-negligable V v can be obtained by scaling the results with 
[r + r v ]. Following this result, the calculations in this paper are 
performed for r = Is -1 and T v = s _1 . Thus, in all the re- 
sults of this paper, the emerging brightness temperature TbAO 
can be replaced with 7\,AQ x [F + T v ] to obtain the results for 
different values of the decay and cross-relaxation rates. 
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Fig. 2. The dependence of the magnetic field strength determined from circular polarization observation, represented by the A F _ F , 
coefficient from Eq.|2l on the magnitude of the velocity gradient AV m for two values of intrinsic thermal line width v t h. Af-f is 
normalized by A° F F ,, which is Af-f for AV m = km s _1 . The solid lines represents masers with different emerging brightness 
temperatures. Note that the noise (< 1%) on the curves is due to the velocity resolution of our model spectra. The results are 
independent of the magnetic field angle 0. 



To solve the equations of state it is assumed that the 
Zeeman frequency shift gQ. is much greater than R, F and 
T v . Then the off-diagonal elements of the density matrix de- 
scribing the molecular states are negligible, greatly simplify- 
ing the calculations. For a magnetic field strength of ~1.0 G, 
gQ. of the dominant hyperfine transitions is ~10 4 s while for 
7bAQ < 10 12 K sr the rate for stimulated emission R < 100 s _1 
(V02). Most of the calculations in this paper are performed for 
a magnetic field strength of B = 1 G. For the highest 7t,AQ, 
the results will be valid for fields down to ~10 mG, when gQ. 
becomes comparable to R. However, the low observed linear 
polarization of the 22 GHz H2O masers (e.g. V02, V05a) indi- 
cates that 7t,AQ is at most 10 11 K sr and typically even lower. 
Thus, the results presented in this paper are valid for B well 
below 1 mG. 

22 GHz H2O masers begin to saturate when R/T > 1 s _1 
and approach full saturation when R/T > 100 s _I . Thus, for our 
models with r = 1 s -1 , saturation starts at 7b AQ > 10 K sr 
and full saturation occurs for T^AD. > 10 12 K sr. The maser 
spectral line profile (through rebroadening) and the polariza- 
tion properties of the maser changes at R = (F + T v ) (e.g. 
NW92). While typically for the 22 GHz H 2 masers this is 
close to when saturation occurs, a large cross-relaxation rate 
means the changes will only occur when the maser has al- 
ready become partly saturated. As T v for the astrophysical 
H2O masers is at most « 5 s for H2O masers at T ~ 1000 K 
and is less than that lower temperatures JAnderson & Watsoi] 
w r is assumed in this paper when discussing unsat- 
urated and saturated masers. It is worth noting that H2O maser 



observations indicate that in actual ity T is between one and tw o 
orders of magnitude less than F v JAnderson & Watsonl fr993). 
However, as our results scale linearly with (r + F v ) this does 
not affect the conclusions in this paper. 

3. Radiative Transfer Modeling Results 

3.1. Circular Polarization 

As shown in NW92 and V02, the V-spectra produced from the 
radiative transfer equations are not anti-symmetric. They are 
also not always directly proportional to the derivative of the 
I-spectrum due to hyperfine interaction as is often assumed. 
However, as was discussed in V02, it is impossible to directly 
observe the intrinsic shape of the V-spectrum due to instrumen- 
tal effects and necessary data calibration steps. Still, the per- 
centage of circular polarization Py can be related to the mag- 
netic field strength using 

^V — (^max ^mirO/Anax 

= 2 -Af-f -^Gaussicose/Aviikms- 1 ]. (2) 

Here V max and V m i n are the minimum and maximum of the V- 
spectrum. 7 max and Avl are the peak flux and the full width 
half-maximum (FWHM) of the I-spectrum respectively. B is 
the magnetic field strength at an angle 8 to the maser propaga- 
tion direction. While in most analyses of maser circular po- 
lariz ation, the coefficient Af -f is taken to be a fixed value 
(e.g iFiebig & Guester]T l989). it was shown in NW92 to de- 
pend on the level of maser saturation and the intrinsic ther- 
mal velocity v t h of the maser. Figure 6 of V02 shows Af-f 
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Fig. 3. a) Derived magnetic field strength B' over the average field strength along the maser (B) as a function of 7bAO. The 
solid and dashed lines correspond to positive and negative magnetic field gradients respectively, b) A' F _ F , as a function of AO for 
different values of emerging brightness temperature 7b AQ. Calculations were performed for Vth = 0.8 km s _1 . 



for different values of v t h as a function of emerging maser 
brightness temperature 7t,AQ. For maser brightness tempera- 
tures r b AQ > 10 9 K sr it was shown in NW92 that the cos# 
dependence of Eq.[2]breaks down effectively also introducing a 
dependen ce on 6 to Af-p . This was later shown in more detail 
in Watson & Wvlcj (2001) for masing involving angular mo- 
mentum / =1-0 and J =2-1 transitions. In V02, figure 7 shows 
the derived magnetic field strength dependence on 6 to Af-f 
for the 22 GHz J =6-5 transition. 

3.1 .1 . velocity gradients 

First, a velocity gradient is introduced by uniformly shifting, in 
Eq.^ the velocity v s of the pump rate Af for each integration 
step along the maser path. By varying the amount of shift, a 
velocity difference AV,„ is created between the start and end of 
the maser amplification path of up to ~1.5 km s _1 . This corre- 
sponds to a velocity gradient of AV„,/S km s _1 m _1 , where S is 
the length of the maser path. As S depends on the exact values 
of the pumping rate Af, while our results, expressed in emerg- 
ing brightness temperature are only dependent on the ratio of 
pumping rates as addressed above, any mention of the velocity 
gradient will henceforth be referring to the value of AV m . In 
Nedoluha & Watson ( 1988) it was shown that the maser splits 
into several distinctly separate narrow features when AV m be- 
comes a few times Vth, each of the features becoming an in- 
dependent maser. The calculations in this paper are therefore 
limited to values of AV m up to ~1.5 km s _1 . 

Fig- n shows the effect of a velocity gradient on the total 
intensity and circular polarization profiles for three different 
intrinsic thermal velocities Vth, produced for the same magnetic 
field strength (1 G). While for low AV,„, the effect of the ve- 
locity gradient on the shape of the I-profile is small, the effect 



on the V-spectrum can still clearly be seen. To examine how 
a velocity gradient influences the magnetic field determination 
from the circular polarization measurements, the normalized 
Ap-f coefficient from Eq.[2]is plotted in Fig.[2]as a function 
of velocity gradient for v t h = 0.6 and 1.5 km s'and emerg- 
ing brightness temperatures T\,AQ. = 10 9 , 10 1() and 10 11 K sr. 
Because of the normalization, the results are independent of 
the magnetic field angle 8. For v t h = 0.6 km s _1 the effect of 
the hyperfine components is clear. For 7bA£2 = 10 9 K sr none 
of the hyperfine components are becoming saturated yet, thus 
Ap-p> is mostly symmetric around AV m = km s _1 . The sym- 
metry breaks down when the hyperfine lines are slowly starting 
to saturate for TbAQ >10 10 K sr, as maser growth is somewhat 
inhibited for velocity gradients in the direction of the weaker 
hyperfine components (F — F' = 6 - 5 and 5 - 4), while a gra- 
dient in the opposite direction allows for stronger maser am- 
plification. As noted above this effect can be postponed further 
into the saturated regime when F v increased. It can be seen in 
the figure, that Af-f typically decreases for increasing |AV,„|, 
implying that in the presence of velocity gradients the mag- 
netic field strength B' derived with Eq. [2] underestimates the 
true B. For small Vth the field strength can sometimes be un- 
derestimated by more than a factor of 2. However, for a nar- 
row range in AV m , Af-f is actually enhanced. In those cases B 
could actually be overestimated by up to ~20%. This effect is 
largest for the more saturated masers. 

3.1 .2. magnetic field gradients 

There are different ways a magnetic field gradient along the 
maser amplification path can occur. The actual field strength 
B can change or the angle between the field and the maser 
propagation direction can vary. Similar to AV m , differences AO 
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Fig. 4. a) Magnetic field angle gradient AO vs. the fractional linear polarization for three different stages of maser saturation and 
6>o = 0°. b) The angle between the maser propagation direction and the magnetic field vs. the fractional linear polar ization for 
different values of emerging maser brightness. The thick solid line denotes the theoretical limit from lGoldreich et alJ dl973h for 
a completely saturated maser. 



and AB are introduced between the start and end of the maser 
path and their influence on the maser polarization is examined. 
Fig. [3^ shows the ratio between the magnetic field B' deter- 
mined with Eq.[2]and the average magnetic field (B) along the 
maser as a function of TbAO when a gradient AB is included. 
The shape of this function does not depend on the actual value 
of AB or on the intrinsic thermal line width v t h. While the maser 
is unsaturated, or more precisely while R < F v , B' is equal to 
(B). However, when the maser saturates, B' will be dominated 
by the average magnetic field in the unsaturated maser region. 

Fig. [3j> shows the change in A' F _ F , versus A9 for different 
TbAQ.. The angle 60, the value of at the start of maser amplifi- 
cation is taken to be 0°. Here A' F _ F , - Af-pcosO, as the actual 
angle cannot be specified due to the inclusion of the gradient 
A6. Comparing the unsaturated maser (TbAQ = 10 9 K sr) with 
the partly saturated maser (I^AQ = 10 11 K sr) it is again appar- 
ent that, as seen in Fig. [3k the unsaturated maser regime pro- 
vides the dominant contribution to the magnetic field strength 
determined with Eq.|3 Comparing Fig.[3p with figure 7 from 
V02 indicates that the angle average (0) in the unsaturated 
maser region determines the observed circular polarization. 

3.2. Linear Polarization 

The linear polarization in the presence of a magnetic field 
was also discussed in NW92. The fractional linear polariza- 
tion is not found to be affected when including a velocity gra- 
dient AV m or magnetic field gradient AB in the calculations. 
Only when including a gradient AO does the linear polarization 
change. In Fig.|4k the fractional linear polarization |<2olM) is 
shown as a function of AO. This can be compared with |<2olM) 
as a function of the angle between the magnetic field and the 



maser propagation direction shown in Fig.|5>. This relation- 
ship, in the limitin g case of a co mpletely saturated maser, was 
solved in lGoldreich et alJ dl973h . In the case of an unsaturated 
maser (TbAQ. < 10 10 K sr), the relationship between \Qq\/Iq 
and the magnetic field angle gradient AO is very similar to the 
relationship between |£2ol/^o and 0, while when the saturation 
level increases the linear polarization is quenched and shifted. 
This indicates that the observed fractional linear polarization 
is determined by at the end of the unsaturated maser regime, 
where the largest amplification occurs, and not by the average 
{0} over the unsaturated maser core. 



4. Conclusions 

Using a maser radiative transfer code which includes the three 
strongest hyperfine components of the 22 GHz H2O 616 - 523 
rotational transition and their magnetic sub-states, the effects 
of velocity and magnetic field gradients on the determination 
of magnetic fields from maser polarization observations have 
been calculated. It was shown that, due to velocity gradients, 
the magnetic field strength determined on H2O masers in CSEs 
(V02, V05a), where v± «0. 8-1.0 km s~'and where velocity 
gradients of ~1.0 km s~'have been observed (V05b), can be 
underestimated by up to a factor of 2. However, for maser 
features with smaller velocity gradients there will be cases 
where B is actually overestimated by ~ 10-20%. The magnetic 
field s determined on the H2O masers in star-form ation regions 
jFiebig & Guestenll989llSarma et alJ200lll2002t) . which have 
higher thermal line widths (v t h ~1.5 km s _1 ) can be underes- 
timated by ~20-40%. These uncertainties can only be over- 
come by including a consistent model for the velocity gradients 
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throughout the H2O maser source when deriving the magnetic 
field strengths derived from circular polarization observations. 

When there is a magnetic field gradients along the maser 
amplification path, B' determined from polarization observa- 
tions of an unsaturated maser (more formally when R < F v ), 
is typically equal to the average field along the maser. When 
the maser is saturated, the B' is dominated by the average mag- 
netic field strength or the average magnetic field angle in the 
unsaturated maser core. The fractional linear polarization on 
the other hand, is not affected by velocity or magnetic field 
gradients. Only a gradient in the angle between the magnetic 
field and the maser propagation axis alters the linear polariza- 
tion. In contrast to the circular polarization, the fractional linear 
polarization is mainly determined by the magnetic field angle 
in the part of the unsaturated maser where the largest ampli- 
fication occurs and not by the angle averaged over the entire 
unsaturated maser region. 
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